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Abstract 

Differential equations which describe the steady state spatial evolution of nonclassical light are 
established using standard quantum field theoretic techniques. A Schrodinger equation for the state 
vector of the optical field is derived using the quantum analog of the slowly varying envelope 
approximation (SVEA). The steady state solutions are those that satisfy the time independent 
Schrodinger equation. The resulting eigenvalue problem then leads to the spatial propagation 
equations. For the degenerate parametric amplifier this method shows that the squeezed state is the 
ground state of the squeezing Hamiltonian. The magnitude and phase of the squeezing parameter 
obey nonlinear differential equations coupled by the amplifier gain constant and phase mismatch. The 
solution to these differential equations is equivalent to one obtained from the classical three wave 
mixing steady state solution to the parametric amplifier with a nondepleted pump. 


1 . Introduction 

The standard approach in quantum optics for dealing with nonclassical light is to introduce a 
normal mode expansion for the field. This naturally leads to time evolution equations for the mode 
amplitudes (creation and annihilation operators) which satisfy the canonical equal time commutation 
relations. Such an approach is well suited for dealing with systems in optical cavities, such as an 
optical parametric oscillator, but is not appropriate for open systems such as a parametric amplifier. 
For the amplifier we are usually interested in the spatial dependence of temporally steady state fields. 
This problem has been treated for quantum fields either by identifying spatial evolution with temporal 
evolution for an interaction time set by the effective interaction length [1-4], or through the use of new 
quantization methods which require nonconventional commutation relations that are generally 
inconsistent with the fundamental equal time commutators for interacting fields [5-9]. In addition, the 
problem has recently been studied using the positive-P representation and stochastic differential 
equations [10]. In this letter we propose an alternative scheme which employs only standard quantum 
field theoretic techniques, and apply this analysis to the degenerate parametric amplifier (DPA) in 
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order to study the spatial evolution of squeezed light. 

In the classical three wave mixing analysis of the DPA, the positive frequency component of the 
electric field is written as a slowly varying time independent envelope multiplying a carrier plane 
wave. Such an ansatz physically represents steady state propagation and leads to spatial differential 
equations for the envelope. Quantum mechanically, this ansatz cannot hold as an operator identity in 
the Heisenberg picture. Upon closer inspection, we see that the requirement of steady state 
propagation is enforced through the choice of state vector for the field. More precisely, steady state 
fields are the stationary states of the appropriate Hamiltonian. The resulting eigenvalue problem is 
then a time independent Schrodinger equation, and leads to spatial differential equations for the 
functions which parameterize the global state vector. We will apply this technique to the DPA in one 
dimension with a classical nondepleted pump. The stationary state solutions are squeezed states 
specified by a set of coupled nonlinear differential equations for the magnitude and phase of the 
squeezing parameter driven by the coupling constant of the three wave mixing. 


2. Steady State Propagation Condition 

Using the classical wave mixing analysis as our guide, we formulate the quantum theory in terms 
of the SVEA. The positive frequency electric field operator is written as an envelope field modulating 
a carrier plane wave of a given polarization E. For systems in which the medium has no transverse 
dependence over the beam diameter and is longitudinally short compared to the Raleigh range, the 
theory becomes essentially one dimensional so that. 


EL +) (x,t) = e 
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( 1 ) 


where A is the beam area, and n(co) is the linear index of refraction. The effective field theory in the 
SVEA formally resembles a nonrelativistic many body theory for a complex scalar field 'F [12], 


mz,D , ^ V,t)] = 8(z-z'). 


( 2 ) 


The Hamiltonian (in the Heisenberg picture) can be written as, 


H — He + H en v +Hi nt 
Hc= flcoj* d 3 x x F t 'F 
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where we have decomposed the "free" field Hamiltonian into the carrier wave Hamiltonian H c , and 
the Hamiltonian governing the free (plus linear) evolution of the envelope H env . The effective 
interaction Hamiltonian for the DPA is given by H^ [14] representing the mixing through the 
nonlinear susceptibility, with pump frequency C0p. The interaction coupling is determined by 


k(z) = ^-g 0 (z)e i< t>( z ), 

0(z) = - + Ak z + ■O(z) 

where go(z) is the standard power gain coupling constant [13] (allowing for possible z-dependence), 
Ak is the phase mismatch at the degenerate frequency, and t3(z) is the remaining phase stemming from 
the pump and X (2) . 

We seek the quantum field state vector corresponding to steady state propagation. Maintaining 
our close analogy with the classical analysis, we remove the carrier wave oscillation from the 
dynamics of this state. Quantum mechanically we use the envelope picture (E.P.) originally 
introduced by Caves [1 1], in which the dynamics of the operators are dictated by the carrier wave 
Hamiltonian, and the states evolve by a Schrodinger equation with the Hamiltonian ^env ^int' ^ 
we denote objects in this picture by the superscript (E), then the dynamical equations in this picture 
are 


(4a) 

(4b) 

(4c) 


= (5a) 

ot 

= <H® + (5b) 

For steady-state applications there will be exact frequency matching between the carrier frequencies of 
the various waves which mix through the nonlinear interaction, so that the Hamiltonian in Eq. (5b) 
will be independent of time. The steady state (ss) solutions are thus identified with the stationary 
solutions to Eq. (5b), 


(H® + H®)KD>® = 


(6) 
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Henceforth, the labels (E) and (ss) will be omitted. Eq. (6) is the key result of our analysis which we 
now apply to the parametric amplifier. 


3. Propagation Solution to the DPA 


It is well known that the output produced during a parametric interaction is a squeezed vacuum 
state [11]. In an oscillator configuration one can restrict attention solely to a single mode of the cavity 
(for a doubly degenerate oscillator). In the amplifier configuration under consideration here, we must 
generalize the squeezing operator to take into account the many modes associated with the envelope 
field. We define a functional squeezing operator, 


S[C) = exp 



(7) 


with z-dependent squeezing parameter £(z) = r(z)e 10 W One can easily show that the Bogoliubov 
transformation generalizes to 


'F(z) = S t [£] x F(z)S[Q = cosh(r(z)) ^(z) - e i0 ( z > sinh(r(z)) 'F + (z). (8) 

The squeezed vacuum. 


i{0> = S[QI0>, 


satisfies the well known eigenvalue equation [15], 


(9) 


5(z)HU> =0, (10a) 

- — V 

¥(z) = S[Q x P(z)S t [Q = cosh(r(z)) T'(z) + e i0 ( z > sinh(r(z)) ^(z). (10b) 

Using the formalism developed in Section 2, we deduce the steady steady solutions by solving 
the time independent Schrodinger equation (6). We take as our ansatz for the eigenvector a squeezed 
vacuum state Eq. (9), i.e. 

(H env + H mt )l{C}> = AJ{Q> . (ii) 

Our goal is to prove Eq. (1 1) and thereby obtain differential equations for the spatial squeezing 
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parameter ^(z). Together, Eqs. (10a) and (11) imply 


[ 


H env + H int ,^(z)Jl(C}>=0 


(12) 


Substituting Eqs. (3) and (10b) into Eq. (12), using Eq. (8), and the relation (4a) for the polar 
decomposition of K gives, 

|jd£ . Ig 0 cos (0-<t>)J + s inh(2r) + go cosh(2r) sin(0-<|>)) Jh^IO) =0. (13) 

Thus, the eigenvalue condition requires the real and imaginary parts of the coefficient of to vanish, 
yielding the desired propagation equations: 


^ = Jrgocos (0-4>) ^ 14a ) 

dz 2 

dfi. = - go coth(2r) sin(0-<f>). (14b) 


The solution to (14) produces a relation for the squeezing parameters (r(z), 0(z)} in terms of the 
experimental parameters { g 0 (z) , <))(z)}. In the next section we give an interpretation of the meaning 
of these equations and analyze their solutions. 


4. Solutions and Interpretation of Propagation Equations 

In solving Eqs. (14) we take as our boundary conditions, 

r(0) = 0 ( 15a ) 

0 ( 0 ) = 0 ( 0 ) ( 15b ) 

to ensure the continuity of the squeezing parameter. In order to get a better understanding of the 
nature of the equations we examine them for various regimes. Consider first the situation in which 
the phase of k, Eq. (4c), is a constant <t>(z) = <|>o (i.e exact phase matching). With boundary 
conditions (15), the solution to (14) is immediately seen to be 


f 


r(z) = I dz' ^go(z'). 


0(z) = <(>o- 


(16a) 

(16b) 
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For a uniform crystal with plane wave pump, this yields the expected result r out = 1-goL (for z>L), 
where L is the length of the crystal. 

In limit of short crystals, Eqs. ( 14) can be expanded for small r to yield 


I ^ = ^gocos (e-<t» | 

I £-§*■<*» f 




(17) 


Thus, in this limit the complex squeezing parameter is just the integral of the interaction coupling 
constant. Henceforth we consider crystals with uniform nonlinearity and a plane wave pump, so that 
g 0 is a constant and (|>(z) = - & + Akz. The solution to Eq. (17) for z>L is of a familiar form, 


Tout — 2^®^ 


sin(AkL/2) 

(AkL/2) 


(18a) 


©out = -| + AkL/2. 


(18b) 


If in addition Ak«l/L the amplitude and phase essentially decouple. However, as the phase 
mismatch increases, the magnitude of the squeezing decreases. Physically, large phase mismatch 
would cause squeezing along wildly different quadratures as the beam propagates through the crystal, 
thereby degrading the net squeezing. 

For crystals containing many gain lengths, we expect strong sensitivity to the phase mismatch 
since the nonlinear coupling between r and 0 causes rapid oscillations that destroy the magnitude of 
the squeezing at some value of Ak. To see this we numerically integrate Eq. (14). In presenting the 
data, we measure all distances in units of the natural length scale l/g 0 . Fig.(l) shows the magnitude 
of the squeezing parameter as a function of z for various phase mismatches. Fig, (2) shows the 
squeezing parameter of downconverted light as function of phase mismatch for various values of the 
crystal length L, plotted simultaneously with the small r expression, Eq. (18a). For short lengths, 

Eq. (18a) is an excellent approximation. However, for crystals many gain lengths long (as can be 
obtained using intense pumps [ 16J) we see a sharp cutoff in the squeezing at the critical value Ak/gQ 
=1. The sharp transition in output when Ak=g 0 is also characteristic of the classical three wave 
mixing solution to the parametric amplifier. This leads us to hypothesize that the solution to Eqs. 
(14) can be obtained from the classical propagation equations. 

The classical solution for the envelope in a DPA is [13] 


A(z) = (p.(z)A(0) + iv(z)A*(0))e iAkz/2 


(19a) 



j_l(z) = cosh(lgz) - i^sinh(lgz) ; v(z) = V 1- \[L(zf = ysinh(^gz) 


(19b) 


g = Vg§- (Ak) 2 . (19c) 

If we replace the c-number fields with creation/annihilation operators A *(())-* at, A(z)->b 
representing "in" and "out" fields respectively, we recover the solution obtained by Caves and Crouch 
[5], and Huttner et. al. [6], 


b = (p(z)a + iv(z)at)e lAkz/2 . 


( 20 ) 


The unitary operator relating a to b, b = U + aU, is a generalized squeezing operator 

U = Rt|l(Arg(|i) + Akz)j s(-isinh _1 (v)), (21) 

where S is the standard squeezing operator and R is the phase rotation operator [11]. This unitary 
operator transforms the vacuum into a squeezed state with squeezing parameter £ = re 10 given by, 

r(z) = sinh 4 (v(z)) 

0(z) = - + Arg(|i(z)) +Akz 

Figs. (1) and (2) are now easily understood in terms of limiting cases of (22a). In fact, direct 
substitution reveals that (22) is indeed the analytic solution to Eq. (14). 

To understand the meaning of this result consider the nature of the interaction more carefully. 
Because we have assumed a classical nondepleted pump, the propagation equations are linear in the 
signal and idler amplitudes . In such a situation one expects the q-number and c-number solutions to 
be equivalent, i.e., the quantum fluctuations should propagate like a classical signal. This provides a 
substantial test on our formalism. The successful use of nonconventional commutation relations by 
previous workers is explained by a combination of good physical intuition and the special linear 
character of this problem. Indeed, Huttner et. al. showed that for this particular interaction the equal-z 
commutation relations are consistent with the fundamental equal-t commutators. However, for a 
general interaction in which the wave mixing equations are truly nonlinear we expect this method to 
fail. For example, if the pump is treated as a quantum mechanical field, then we find that the equal-z 
commutation relations are no longer consistent with the fundamental equal-t commutators. Thus, for 
more complex systems there will be qualitative differences between the classical and quantum 
propagation solutions. The formalism we have presented here is sufficiently general to handle these 
situations. 


(22a) 

(22b) 
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le as a function of propagation distance z measured in gain lengths g 0 , 
es of phase mismatch Ak: (a) Ak=0.1g 0 , (b) Ak=g 0 , (c) Ak=2g 0 . 
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BgJL Squeezing at the output face as a function of Ak: (a) L=2g 0 ,(b) L=10g 0 , (c) L=50g 0 . The 
solid curve is the solution to Eq. (14), the dotted curve is given by the asymptotic result Eq. (18a). 
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